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ABSTRACT 

We describe the results of an extremely deep, 0.28 deg^ survey for z = 3.1 
Lya emission-line galaxies in the Extended Chandra Deep Field South. By using 
a narrow-band 5000 A filter and complementary broadband photometry from 
the MUSYC survey, we identify a statistically complete sample of 162 galaxies 
with monochromatic fluxes brighter than 1.5 x 10""^^ ergs cm~^ s~^ and observers 
frame equivalent widths greater than 80 A. We show that the equivalent width 
distribution of these objects follows an exponential with a rest-frame scale length 
of wq = 76]^g^ A. In addition, we show that in the emission line, the luminosity 
function of Lya galaxies has a faint-end power-law slope of a = — 1.49]']q 34, a 
bright-end cutoff of logL* = 42.641°;^^, and a space density above our detection 
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thresholds of 1.46 ± 0.12 x 10~^ hjQ^ galaxies Mpc~^. Finally, by comparing the 
emission-line and continuum properties of the LAEs, we show that the star- 
formation rates derived from Lya are ~ 3 times lower than those inferred from 
the rest-frame UV continuum. We use this offset to deduce the existence of a 
small amount of internal extinction within the host galaxies. This extinction, 
coupled with the lack of extremely-high equivalent width emitters, argues that 
these galaxies are not primordial Pop III objects, though they are young and 
relatively chemically unevolved. 

Subject headings: cosmology: observations - galaxies: formation - galaxies: high- 
redshift - galaxies: luminosity function 



Introduction 



The past decade has seen an explosion in our ability to dete ct and study ^ > 3 galax- 
ies and probe the history of star formation in the universe (e.g., iMadau et al.lll996l ). This 
has been mostly due to the development of the Lyman-break technique, whereby high red- 
shift galaxies are id entified via a flux discontinuity caused by Lyman-limit absorption (see 
Steidel et al.lll996al Jbl). By taking deep broadband images, and searching for U, B, and V- 
band dropouts, astronomers have been able to explore large-scale stru cture and determ ine 
the properties of bright (L > 0.3L*) galaxies between z ~ 3 and z ~ 5 (iGiavaliscol |2002| ) . 



The stunning success of the Lyman-break technique stands in contrast to the initial 
results of Lyg emission-line observations. The failure of the first gener ation of these sur- 



veys (e.g., 



De Propis et al.lll993l : [Thompson. Diorgovski. fc 



internal extinction in the target galaxies (IMeier fc Terlevichlll98ll ). Since Lya photons are 



raugerlll995l ) was attributed to 



resonantly scattered by interstellar hydrogen, even a small amount of dust can reduce the 
emergent emission-line flux by several orders of magnitude. 

Fortunately, Ly a surveys have rec ently undergone a resurgence . ^Starting with the 
Keck observations of ICowie fc Hul (Il998l ) and lHu. Cowie. fc McMahonI (Il998l ). narrow-band 
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t al 2001). z ~ 3.1 fCiarduUo et al. 2002: 


Havashino et al. 


2004: 


Venemans et al. 2005: 


Gawiser et al. 


2006a). ^ ~ 3.7 f 


Fuiita et a . 


200 


3), z ~ 4.5 


(Rhoads et al. 2000 


). ^ ~ 4.9 (Ouchi et al. 


2003h. z 5.7 (Rhoads et al. 


2OO3I; 


Ajiki et al. 


2OO3I; 


Tapken et al. 


2006). and z ~ 6.5 (Kodaira et al. 20031: 


Taniguchi et al. 


2005h. The 



discovery of these high-redshift Lya emitters (LAEs) has opened up a new frontier in as- 
tronomy. At 2; > 4, LAEs are as easy to detect than Lyman-break galaxies (LBG), and, by 
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z > Q, they are the only galaxies observable from the ground. By selecting galaxies via their 
Lya emission, it is therefore possible to probe much further down the galaxy continuum 
luminosity function than with the Lyman-break technique, and perhaps identify the most 
dust-free object s in the universe. In addition, by using Lya emitters as tracers of large- 
scale structure (ISteidel et al.ll2000l : IShimasaku et al.l 120041 ). it is possible to efficiently probe 
the expansion history of t he universe with a minimum of cosrnologi cal assumptions (e.g., 
Blake fc GlazebrooklEoOsI : ISeo fc Eisensteinlboosl : iKoehler et aPboOTh . 



Here, we describe the results of a deep survey for Lya emission-line galaxies in a 
0.28 deg^ region centered on the Extended Chandra Deep Field South (ECDF-S). This 
region has an extraordinary amount of complementary data, including high-resolution op- 
tical images from the Hubble Space Te lescope via the Great Observatories Origins Deep 



Survey (GOODS; iGiavalisco et a 



20041) and the Galaxy Evolution from Morphology and 



SEDs program (GEMS mix et al.ll2004l ). deep groundb ased U BVRIzJHK p hotometry from 
the Multiwavelength Survey by Yale-Chile (MUSYC; iGawiser et al.ll2006bl ). mid- and far- 
IR observations froni Spitzer, GOODS and MUSYC, and deep X-ray data from Chandra 
(iGiacconi et al.l 120021 : [Alexander et al.ll2003l : iLehmer et al.ll2005l ). In Section 2, we describe 
our observations, which include over 28 hours worth of exposures through a narrow-band filter 
on the CTIO 4-m telescope. We also review the techniques used to detect the emission-line 
galaxies, and discuss the difficulties associated with analyzing samples of LAEs discovered 
via fast-beam instruments. In Section 3, we describe the continuum properties of our Lya 
emitters, including their rest-frame mioso — ""^1570 colors, and compare their space density 
to that of Lyman-break galaxies. In Section 4, we examine the LAE's equivalent width 
distribution and show t hat our sampl e cont ains very few of the extremely-high equivalent 
width objects found by iDawson et al.l (120041 ) at z = 4.5. In S ection 5, we pres ent the Lya 
emission-line luminosity function, and give values for its best-fit ISchechterl (119761 ) parameters 
and normalization. In Section 6, we translate these Lja fiuxes into star-formation rates, and 
consider the properties of LAEs in the context of the star-formation rate (SFR) history of 
the universe. We conclude by discussing the implications our observations have for surveys 
aimed at determining cosmic evolution. 



M 



For our analysis, we adopt a ACDM cosmology with Hq = 70 km s ^ Mpc ^ (/170 = 1), 
= 0.3, and = 0.7. At z = 3.1, this implies a physical scale of 7.6 kpc per arcsecond. 



2. Observations and Reductions 



Narrow-band observations of the ECDF-S were performed with the MOSAIC II CCD 
camera on the CTIO Blanco 4-m telescope. These data consisted of a series of 111 exposures 
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taken over 16 nights through a 50 A wide full-width-half-maximum (FWHM) A5000 filter 
(see Figured]). The total exposure time for these images was 28.17 hr; when the effects of 
dithering to cover for a dead CCD during some of the observations are included, the net 
exposure time becomes ~ 24 hr. The total area covered in our survey is 998 arcmin^; after 
the regions around bright stars are excluded, this area shrinks 993 arcmin^. The overall 
seeing on the images is I'.'O. A log of our narrow-band exposures appears in Table [H 

The procedures used to reduce the d ata, identify line emitte rs, a nd measure their bright- 
nesses were identical to those detailed in lCiarduUo et al.l (120021 ) and iFeldmeier et al.l (120031 ). 
After de-biasing, flat-fielding, and aligning the data, our narrow-band frames were co-added 
to create a master image that was clipped of cosmic rays. This frame was then comp ared to 
a deep B+V continuum image provided by the MUSYC survey (iGawiser et al.ll2006bl ) in two 
different ways. First, the DAOFIND task within IRAF was run on the summed narrow-band 
and continuum image using a series of three convolution kernels, ranging from one match- 
ing the image point- spread-function (PSF), to one ~ 3 times larger. This created a source 
catalog of all objects in our field. These targets were then photometrically measured with 
DAOPHOT's PHOT routine, and sources with on-band minus continuum colors less than 
—1.03 in the AB system were flagged as possible emission-line sources (see Figure [2]). At 
the same time, candidate LAEs were also identified by searching for positive residuals on a 
"difference" image made by subtracting a scaled version of the B+V continuum image from 
the narrow-band frame. In this case, the DAOFIND algorithm was set to flag all objects 
brighter than fou r times the local stand ard deviation of the background sky (see Figure [3]) • 
As pointed out by IFeldmeier et al.l (120031 ) . these two techniques complement each other, since 
each detects objects that the other does not. Specifically, < 10% of galaxies were missed 
by the color-magnitude method due to image blending and confusion, but found with the 
difference method. Conversely, objects at the frame limit that were lost amidst the increased 
noise of the difference frame, could still be identified via their on-band minus off-band colors. 



Finally, because we intentionally biased our DAOFIND parameters to identify faint 
sources at the expense of false detections, each emission-line candidate was visually inspected 
on the narrow-band, B + V continuum, and difference frames, as well as two frames made 
from subsamples of half the on-band exposures. This last step excluded many false detections 
at the frame limit, and left us with a sample of 259 candidate LAEs for analysis. 

Once found, the equatorial positions of the candidate emission-line gala xies were de- 



rived with respect to the reference stars of the USNO-A 2.0 astrometric catalog (IMonet et al. 
19981 ). The measured residuals of the plate solution were ~ 0'.'2, a number slightly less than 



the 0'.'25 external error associated with the catalog. Relative narrow-band magnitudes for 
the objects were derived by first measuring the sources with respect to field stars using an 
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aperture slightly greater than the frame PSF. Since most of the galaxies detected in this 
survey are, at best, marginally resolved on our 1" images, this procedure was sufficiently 
accurate for our purposes. We then obtained standard AB magnitudes by comparing large 
aperture photometry of the field sta rs to similar measurements of the spectrophotometric 
standards Feige 56 and Hiltner 600 (IStond 119771 ) taken on three separate nights. The dis- 
persion in the photometric zero point computed from our standard star measurements was 
0.03 mag. 



2.1. Derivation of Monochromatic Fluxes 



The fast optics of wide-field instruments, such as the MOSAIC camera at the CTIO 4-m 
telescope, present an especially difficult challenge for narrow-band imaging. The transmission 
of an interference filter depends critically on the angle at which it is illuminated: light 
entering at the normal will const ructively/ destructively int erfere at a different wavelength 
than light coming in at an angle (lEather fc Reasoned 1 19691 ). As a result, when placed in a 
fast converging beam, an interference filter will have its bandpass broadened and its peak 
transmission decreased by a substantial amount. This effect is important, for without precise 
knowledge of the filter bandpass, it is impossible to derive accurate monochromatic fiuxes 
or estimate equivalent widths. 

To derive the filter transmission, we began with the throughput information provided by 
the CTIO observator}0. This curve, which represents the expected transmission of the [O III] 
interference filter in the f/3.2 beam of the Blanco telescope, was computed by combining 
laboratory measurements of the filter tippe d at several differen t angles from the incoming 
beam (for a discussion of this procedure, see lJacoby et al.lll989l ). We then shifted this curve 
2 A to the blue, to compensate for the thermal contraction of the glass at the telescope, and 
compared this model b andpass to the measured emission-line wavelengths obtained from 
follow-up spectroscopy (ILira et al.ll2007l ). Interestingly, redshift measurements of 72 galaxies 
detected in three independent MUSYC fields confirm the shape of the filter's transmission 
curve, but not its central wavelength: accordin g to the spectroscop y, the mean wavelength 
of the filter is 10 A bluer than given by CTIO ( Gawiser et al.l 120071 ). Examining the source 
of this discrepancy is beyond the scope of this paper. However, the data do confirm that, 
when placed in the beam of the CTIO 4-m prime focus MOSAIC camera, the bandpass of 
the CTIO [O III] interference filter is nearly Gaussian in shape. This bandpass is reproduced 
in the left-hand panel of Figure |H 



^http://www. ctio.noao.edu/instruments/FILTERS/indcx. html 
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This non-square bandpass has important consequences for the analysis of large samples 
of emission-line galaxies. The first of these involves the definition of survey volume. Because 
the transmission of the filter declines away from the bandpass center, the volume of space 
sampled by our observations is a strong function of line strength. This is illustrated in the 
center panel of Figure |H Objects with bright line emission can be detected even if their 
redshifts place Lya in the wings of the filter, hence the volume covered for these objects is 
realtively large. Conversely, weak Lya sources must have their line emission near the center 
of the bandpass to be observable. As a result, the "effective" volume for our integrated 
sample of galaxies is a function of the galaxy emission-line luminosity function. 

A second concern deals with the sample's flux calibration. In order to compare the flux 
of an emission-line object to that of a spectrophotometric standard star (i.e., a continuum 
source) one needs to know both the filter' s integral transmission a nd its monochromatic 



transmission at the wavelength of interest (jJacoby et al.l 119871 . Il989l ) . When observing ob 



jects at known redshift, the latter requirement is not an issue. However, when measuring 
a set of galaxies which can fall anywhere within a Gaussian-shaped transmission curve, the 
transformation between an objects' (bandpass-dependent) AB magnitude and its monochro- 
matic flux is not unique. In fact, if we assume that galaxies are (on average) distributed 
uniformly in redshift space, then the number of emission-line objects present at a given 
transmission, T, is simply proportional to the amount of wavelength associated with that 
transmission value. Consequently, the observed distribution of emission-line fluxes will be 
related to the true distribution via a convolution, whose (unity normalized) kernel, G{T), is 



G{T)dT 



dT 



dry + 

blue 



dX 



dT 



dT 



red 



where the first term describes the filter's response blueward of the transmission peak and the 
second term gives the response redward of the peak. The center panel of Figure H] displays 
this kernel for the filter used in our survey. The curve shows that for roughly half of the 
detectable galaxies in our field, the effect of our filter's non-square bandpass is minimal. 
However, for the other ~ 50% of galaxies, the shape of the bandpass is extremely important, 
and the inferred fluxes for some objects can be off by over a magnitude. 

Any analysis of the ensemble properties of our LAEs must consider the full effect that 
the non-square bandpass and the odd-shaped convolution kernal has on the sample. We do 
this in Sections 4 and 5. However, one often wants to quote the monochromatic flux and 
equivalent width for an individual Lja emitter. To do this, we need to adopt an appropriate 
"mean" value for the transmission of our filter. The most straightforward way to define this 
number is via the filter's peak transmission. This is where the survey depth is greatest, and 
choosing T^ax is equivalent to assigning each galaxy its "most probable" monochromatic flux. 
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Unfortunately, by defining the transmission in this way, we underestimate the flux from all 
galaxies whose line emission does not fall exactly on this peak. Alternatively, we can attempt 
to choose a transmission which globally minimizes the flux errors of all the galaxies detected 
in the survey. This can be done by weighting each transmission by the number of galaxies one 
expects to observe at that wavelength: the greater the transmission, the deeper the survey, 
and the more galaxies present in the sample. The difficulty with this "expectation value" 
approach is that it requires prior knowledge of the distribution of emission-line fluxes, which 
is one of the quantities we are attempting to measure. That leads us to a third possibility: 
to approximate the filter's expectation value using some "characteristic" transmission, Tc, 
which is independent of the galaxy luminosity function, but still takes the filter's changing 
transmission into account. The arrow in Figure H] identifies the transmission we selected as 
being characteristic of the filter; the justification for this value is presented in Section 5. We 
emphasize that Tc is only a convenient mean that enables us to quote the likely emission- 
line strengths of individual galaxies. When analyzing the global properties of an ensemble 
of LAEs, the full non-Gaussian nature of the filter's convolution kernel must be taken into 
account. 

Using this transmission and our knowledge of the filter curve, we converted the galaxies' 
AB magnitudes to monochromatic fiuxes at A = 5000 A via 

F5000 = 3.63 X 10-^° io-AB/2.5 . c [T^ 



where -F5000 is given in ergs cm ^ s ^ (jjacoby et al.l 119871 ). Equivalent widths then followed 



via 

EW = 1^ - AA (3) 

JB+V 

where Jb+v is the objects' AB fiux density in the B + V continuum image, and AA, the 
FWHM of the narrow-band filter, represents the contribution of the galaxy's underlying 
continuum within the bandpass. Both these equations are only applicable to objects whose 
line emission dominates the continuum within the narrow-band filter's bandpass. Since we 
are limiting our discussion to galaxies with narrow-band minus broad-band AB magnitudes 
more negative than —1.03, this approximation is certainly valid. However, we do note that 
by using Tc instead of Tmaxj "we are intentionally overestimating the fiux and equivalent 
width of some galaxies, in order to minimize the errors in others. So, while the application 
of Tc formally translates our Am = —1.03 criterion into a minimum emission-line equivalent 
width of 90 A, galaxies with emission-lines that fall near the peak of the filter transmission 
function can have equivalent widths that are ~ 12% smaller. This implies that the absolute 
minimum equivalent width limit for our sample of LAEs is 80 A. 
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2.2. Sample of LAE Candidates 



Tables |2] and [3] give the coordinates of each candidate emission-line galaxy, along with 
its inferred monochromatic flux and equivalent width. In total, 259 objects are listed, 
though many are beyo nd the limit of our com pleteness. To determine this limit, we fol- 
lowed the procedures of iFeldmeier et al.l (120031 ) and added 1,000,000 artificial stars (2000 at 
a time) to our narrow-band frame. By re-running our detection algorithms on these modi- 
fied frames, we were able to compute the fiux level below which the object recovery fraction 
dropped below the 90% threshold. This value, which corresponds to a monochromatic fiux 
of 1.5 X 10~^^ ergs cm~^ s~^ (logFsooo = —16.82) is our limiting magnitude for statistical 
completeness; 162 galaxies satisfy this criterion. 

Before proceeding further with our analysis, we performed one additional check on 
our data. To eliminate obvious AGN from our sample, we cross-correlated our catalog of 
emission-line objects with t he lists of X-ray sourc es found in the 1 Msec exposure of the 



Chandra Deep Field South ( Alexander et al. 2003h. and the four 250 ksec ex posure of the 



Extended Chandra Deep Field South (ILehmer et al.l l2005l : IVirani et al.l 120061 ) . Two of our 
LAE candidates were detected in the X-ray band. The first, which is our brightest Lya 
emitter, has a 0.5 - 8 keV fiux of 3.4 x 10"^^ ergs cm~^ s ~^ ('i.e., Lx ' y 2.8 x 10'^'' /i7o~^ ergs s"^ 
ai z = 3.1) and exhibits C IV emission at 1550 A (ILira et al.l 120071 ). The other is an 
interloper: a z = 1.6 AGN detected via its strong C III] line at 1909 A. For the remaining 
160 objects that were not detected individually in the X-ray band, we used stacking analyses 
to constrain their mean X-ray power output (see iLehmer et al.l 120071 . for details). We find 
that the stacked X-ray signal, which corresponds to a ~ 40 ksec effective exposure on an 
average LAE, does not yield a 3a detection in any of three X-ray bandpasses (0.5-8.0 keV, 
0.5-2.0 keV, and 2-8 keV). These results imply a 3a upper-limit of ~ 3.8 x 10^^ h^Q ergs s~^ 
on the mean 0.5-2.0 keV luminosity for our LAEs, which demonstrates that f ew of our Lyg 



sources harbor low-luminosity AGN. Similarly, if we use the conversion of iRanalli et al. 



(120031 ). we can translate this X-ray non-detection into an upper-limit for a typical LAE's 
star- formation rate. This limit, 85 hjQ~'^ Mq yr~^, is roughly an order of magnitude greater 
than the rates inferred from the objects' Lya emission or UV continua (see Section 6). 

For the remainder of this paper, we will treat our z = 3.1 X-ray source as AGN and 
exclude it from the analysis. This leaves us with a sample of 160 objects, which we assume 
are all star-forming galaxies. We note that, because all of our objects have equivalent widths 
greater than 80 A, they are unlikely to be [O II] emitters. At 2; ~ 0.34, our survey volume 
is only ~ 730 /i 70"^ Mpc ^, which, through th e luminosity functions of iHogg et al.l (119981 ) . 
Gallego et al.l (120021 ). and iTeplitz et al.l (120031 ) . implies a total population of between ~ 20 



and ~ 200 [O II] emission-line galaxies above our completeness limit. Since less than 
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of th ese objects will have rest frame equivalent widths greater than ~ 60 A ( Hogg et al. 



19981 ). the number of [O II] interlopers in our sample should be negligible. This estimate 



is confirmed by follow-up spectroscopy: of the 52 LAE candidates observed with sufficient 



signal-to-noise for a re dshift determination, all are confirmed Lya emitters (iGawiser et al. 



2006al : lLira et all 120071 ). 



Figure [5] shows the spatial distribution of the LAEs above our completeness limit. The 
sources are obviously clustered, falling along what appear to be "walls" or "filaments" . The 
GOODS region has a below-average number of 2; = 3.1 Lya emitters, and there are almost 
no objects in the northwestern part of the field. Gonversely, the density of LAEs east and 
northeast of the field center is quite high. This type of data can be an extremely powerful 
probe of cosmologica l history, but we will defer a discussion of this topic to a future paper 
dOawiser et"aDl2007h . 



3. The Continuum Properties of the Emitters 



To investigate the continuum properties of our Lya emitters, we r neasured the brightnes s 
of each LAE on the broadband UBVR images of the MUSYG survey (IGawiser et al.ll2006bl ). 
Since the catalog associated with this dataset has a 5 a detection threshold of U = 26.0, 
B = 26.9, V = 26.4, and R = 26.4, our knowledge of the LAEs' positions (obtained from 
the narrow-band frames) allows us to perform photometry well past this limit. Figure [6] 
displays the B — R color-magnitude diagram for 88 of the LAEs brighter than Rab = 27.25. 
The diagram, which shows the galaxies' rest-frame continua at 1060 and 1570 A, has several 
features of note. 

The first involves the color distribution of our objects. According to the figure, LAEs 
with i?-band magnitudes brighter than R = 25 have a median color oi B — R = 0.53. 



This value agrees with the blue colors found by IVenemans et al.l ( 120051 ) for a sample of Lya 
emitters at z = 3.13, and is the value expected for a ~ 10^ yr old stellar system evolving with 
a constant star- formation rate ( Fuiita et al. l20o"3l: Bruzual fc Gharlot 2003). This median 



color is also consistent with the results of lGawiser et al.l (l2006al ). who stacked the broad-band 
fluxes of 18 spectroscopicaUy confirmed z = 3.1 LAEs and showed that the typical age of 
these systems is between 0.01 < t < 2 Gyr. It does, however, stand in marked contrast to 
the results of IStiavelli et al.l (120011 ). who claimed that Lja emitters at z = 2.4 are very red 
(5 — / ~ 1.8). The blue colors of our galaxies confirm their nature as young, star- forming 
systems. There is no evidence for excessive reddening in these objects, and if the galaxies 
do possess an underlying population of older stars, the component must be quite small. 
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On the other hand, as the LAE color distribution indicates, Lya emitters are not, as 
a class, homogeneous. At R = 25, the MUSYC B — R colors have a typical photometric 
uncertainty of (Tb~r = 0.25 mag. This contrasts with the observed color dispersion for our 
galaxies, which is ~ 0.4 mag for objects with R < 25. Thus, there is at least a ~ 0.3 mag 
scatter in the intrinsic colors of these objects. Either there is some variation in the star- 
formation history of Lya emitters, or dust is having an effect on the emergent colors. 

Finally, it is worth emphasizing that our Lya emitters are substantially fainter in the 
continuum than objects found by the Lyman-break techn ique. At ^ ~ 3, L* galaxies have 
an apparent magnitude of i? ~ 24.5 (ISteidel et al.l Il999l ) a nd ground-based Lyman-break 
surveys typically extend only ~ 1 mag beyond this value (see lGiavaliscdl2002l . for a review). 
Furthermore, spectroscopic surveys of LBG candidates rarely target galaxies fainter than 
R = 24. In our emission-line sample, the median continuum magnitude is R ^ 26.7, and 
many of the galaxies have aperture magnitudes significantly fainter than R ~ 28. In general, 
LAEs do inhabit the same location as LBGs in the U-V vs. V-R color-color space (see 
Figure [7j), but their extremely faint continuum sets them apart. 

This is also illustrated in Figure [S], which compares the rest frame 1570 A luminosity 
function of our complete sample of Lya emitters (those with monochromatic fluxes greater 
than 1.5 x 10"^'' ergs c m~^ s~^) with the rest-frame 1700 A luminosity function of 2; = 3.1 
Lyman-break galaxies (ISteidel et al.lll999l ). When plotted in this way, our sample of LAEs 
appears incomplete, since for R > 26.5, only the brightest emission-line sources will make it 
into our catalog. The plot also implies that at z = 3.1, R < 25.5, Lya emitters are ~ 3 times 
rarer than comparabl y bright Lyrnan-bre ak galaxies. Since this ratio is virtually identical 
to that measured by ISteidel et al.l (I2OOOI ) within an extremely rich z = 3.09 protocluster, 
this suggests that the number is not a strong function of galactic environment. But, most 
strikingly, our observations demonstrate the Lya emitters sample the entire range of the 
(UV-continuum) luminosity function. The median UV luminosity of LAEs in our sample 
is < 0.2L*, and the faintest galaxy in the group is no brighter than ~ 0.02L*. Just as 
broadband observations detect all objects at the bright-end of the continuum luminosity 
function, but sample the entire range of emission-line strengths, our narrow-band survey 
finds all the brightest emission-line objects, but draws from the entire range of continuum 
brightness. 



4. The Lya Equivalent Width Distribution 



Before examining the emission-line properties of our dataset, we need to correct for 
the observational biases and selection effects that are present in the sample. Since the 
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data were taken in a fast-beam through a filter with a non-square bandpass, these effects 
are substantial. Continuum measurements, of course, are unaffected by the peculiarities 
of a narrow-band filter, but the distribution of monochromatic fluxes can be significantly 
distorted. Specifically, the observed flux distribution will be the convolution of the true 
distribution with the following two kernels: 

The Photometric Error Function: The random errors associated with our narrow-band 
photometry vary considerably, ranging from ~ 0.02 mag at the bright end, to ~ 0.2 mag near 
the completeness limit (see Table Hj). These errors will scatter objects from heavily populated 
magnitude bins into bins with fewer objects, and flatten the slope of the luminosit y function. 



Beca.use th e change in slope goes as the square of the measurement uncertainty (lEddington 



19131 . Il940l ). the effect of this convolution is most important for objects near the survey limit. 



The Filter Transmission Function. As described in Section 2.1, the narrow-band filter 
used for this survey has a transmission function that is nearly Gaussian in shape. This 
creates an odd-shaped convolution kernel (the right panel of Figure H]), which systematically 
decreases the measured line-emission of objects falling away from the peak of the transmission 
curve. Moreover, because the objects' equivalent widths are also reduced by this bandpass 
effect, some fraction of the LAE population will be lost from our EW > 80 A sample. 
The result is that the normalization of this filter transmission kernel is not unity. Instead, 
it depends on the intrinsic equivalent width distribution of the galaxies, since that is the 
function that defines the fraction of galaxies (at each redshift) which can still make it into 
our sample. 

These effects are illustrated in the top panel of Figure [9l which displays a histogram of 
the rest-frame equivalent widths for our candidate Lya galaxies. As the dotted line shows, the 
data appear to be well fit by an exponential that has an e-folding length of Wobs = 214^^5 A. 
However, because the bandpass of our narrow-band filter is more Gaussian-shaped than 
square, the line-strengths of many of the galaxies have been systematically underestimated. 
In fact, the true distribution of equivalent widths is broader than that measured: when we 
perform a maximum-likelihood analysis using a series of exponential laws, convolved with 
the filter bandpass and photometric error kernels, we obtain a most-likely scale length of 
Wohs = 311^33 A, or wq = 76^g^ A in the rest frame of the sample. 



Such a distribution is quite different from that reported by iMalhotra fc Rhoadd (120021 ). 
In their survey of 150 z = 4.5 Lya emitters, ~ 60% of the objects had extremely high rest- 
frame equivalent widt hs, i.e., EWq > 240 A. Since stellar populat ion models, such as those by 



Chariot fc Falll (119931 ) cannot produce such strong hne-emission, iMalhotra fc Rhoadd (120021 ) 
postulated the presence of a top-heavy initial mass function and perhaps the existence of 
Population III stars. However in our sample, only 3 out of 160 LAEs (~ 2%) have observed 
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rest-frame equivalent widths greater than this 240 A hmit. Even when we correct for the 
effects of our filter's non-square bandpass, the fraction of stro ng line-emitter s does not exceed 
~ 12%. This is less than the ~ 20% value estimate d by iDawson et al.l (12004 ) via Keck 
spectroscopy of a subset of iMalhotra fc Rhoadsl (120021 ) objects. Thus, at least at z ~ 3.1, 
there is no need to invoke a skewed initial mass function to explain the majority of our LAEs. 

The equivale n t widt h distribution of Figure [9] also differs dramatically from that found 
by IShapley et al.l (120031 ) for a sample of 2; ~ 3 Lyman-break galaxies. In their dataset, 
rest-frame equivalent widths e-fold with a scale-length of ~ 25 A, rather than the ~ 75 A 
value derived from our LAE survey. This difference is not surprising given that the former 
dataset is selected to be bright in th e continuum, while the latter is chosen to be strong 
in the emission-line. Moreover, when IShapley et al.l (120031 ) analyzed the ~ 25% of Lyman- 
break galaxies with rest-frame equivalent widths greater than 20 A, they found a correlation 
between line strength and continuum (i?-band) magnitude, in the sense that fainter galaxies 
had higher equivalent widths. We see that same trend in our data, but it is largely the result 
of a selection effect. (Faint galaxies with low equivalent widths fall below our monochromatic 
flux limit.) A comparison of emission-line flux with equivalent width for our statistically 
complete sample shows no such correlation. 

The lower two panels of Figure [9] demonstrate this another way. In the diagram, our 
sample of LAEs is divided in half, with the middle panel showing the equivalent width 
distribution for objects with monochromatic Lya luminosities greater 2 x 10^^ ^7o~^ ergs s~^, 
and the bottom panel displaying the same distribution for less luminous objects. As the 
figure illustrates, the distribution of equivalent widths is relatively insensitive to the absolute 
brightness of the galaxy. To first order this is expected, since both the UV continuum and 
the Lya emission-line flux are driven by star formation. However, one could imagine a 
scenario wherein the amount, composition, and/or distribution of dust within the brighter 
(presumably more-metal rich) Lya emitters differs from that within their lower- luminosity 
counterparts. Since the effect of this dust on resonantly-scattered Lya photons is likely 
to be different from that on continuum photons, this change in extinction can theoretically 
produce a systematic shift in the distribution of Lya equivalent widths. There is no evidence 
for such a shift in our data; this constancy argues against the importance of dust in these 
objects. 



5. The Lya Emission-Line Luminosity Function 



Figure [TO] shows the distribution of monochromatic fluxes for our sample of emission- 
line galaxies. The function looks much like a power law, with a faint-end slope of a ~ —1.5 
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that steepens as one moves to brighter luminosities. However, to quantify this behavior, 
we once again have to correct the observed flux distribution for the distortions caused by 
photometric errors and the non-square bandpass of the filter. In addition, we must also 
consider the censoring effect our equivalent width cutoff has on the data: some line emitters 
whose redshifts are not at the peak of the filter transmission function will fall out of our 
sample completely. 

To dea l with these effects, we fit the observed distribution of Lya er nission-line fluxes to 
Schechter Jl976h fu nction via the method of maximum likelihood (e.g.. lHanes fc Whittakerl 



1987 



CiarduUo et al. I ll989h . We applied our two convolution kernels (including the equivalent 
width censorship) to a series of functions of the form 

(t){L)d{L/L*) oc {L/L*fe-^/^'d{L/L*) (4) 

treated each curve as a probability distribution (i.e., with a unity normalization), and com- 
puted the likelihood that the observed sample of Lya fluxes is drawn from the resultant 
distribution. The results for the three parameters of this fit, a, logL*, and A^, the integral 
of the Schechter function down to our limiting flux (in units of galaxies Mpc~^), are shown 
in Figure [TTl Table O lists the best-fitting parameters, along with their marginalized most- 
likely values and uncertainties. For completeness. Table [5] also gives the value of 0* which is 
inferred from our most likely solution. As expected, the plots illustrate the familiar degen- 
eracy between L* and a: our best-fit solution has a ~ —1.5, but if L* is forced to brighter 
luminosities, a decreases. The contours also demonstrate an asymmetry in the solutions, 
whereby extremely bright values of L* are included within the 3 a contours of probability, 
but faint values of the same quantity are not. 

But perhaps the most interesting feature of the analysis concerns the effective volume 
of our survey. As in Section 2.1, the amount of space sampled by the observations depends 
critically on each galaxy's Lya luminosity and equivalent width. Bright line-emitters with 
large equivalent widths can be identified well onto the wings of the filter, hence the survey 
volume associated with these objects is relatively large. Conversely, weak line-emitters, and 
objects with small equivalent widths can only be detected if they lie at the peak of the filter 
transmission curve. Thus, the survey volume for these objects is quite small. The effective 
volume for our observations is therefore a weighted average, which depends on the intrinsic 
properties of entire LAE sample. 

This average can be computed from the data displayed in Figure [TTl According to the 
figure, the space density of galaxies with emission-line brighter than 1.5 x 10"^"^ ergs cm~^ s~^ 
(i.e., 1.3x10''^ /iTo"^ ergs s~^) is extremely well-defined, 1.46±0.12xl0~^ /lyo^ galaxies Mpc"^. 
Since this measurement comes from the detection of 160 galaxies brighter than the complete- 
ness limit, the data imply an effective survey volume of ~ 1.1 x 10^ ^7o~^ Mpc'^. This is not 
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the volume one would infer from the interference filter's full-width at half-maximum: it is 
25% smaller, or roughly the full-width of the filter at two-thirds maximum. 

This difference is illustrated in Figure [101 The points show the space density of Lya 
galaxies one would derive sirnply by using the filter's FWHM to define the survey volume; the 



solid line gives the ISchechterl (119761 ) function which best fits the data. The offset between the 
solid line and the dashed line, which represents the function after the application of the two 
convolution kernels, confirms the need for careful analysis when working with narrow-band 
data taken through a non-square bandpass. 

The results of our maximum-likelihood calculation also suggest a simple definition for 
the effective transmission for our filter. As described in Section 2.1, a "characteristic" trans- 
mission is needed to convert the (bandpass-dependent) AB magnitude of an individual galaxy 
to monochromatic Lya fiux. Rather than use the maximum transmission (which would un- 
derestimate the fiux of all galaxies not at the filter peak), or adopt some complicated scheme 
which involves iterating on the luminosity function, one can simply choose the filter's mean 
transmission within some limited wavelength range. Based on the results above, the filter's 
full- width at two-thirds maximum seems an appropriate limit. This transmission, which is 
indicated by the arrow in Figure HI is the value used to derive the fiuxes and equivalent 
widths of Tables 121 and [3l If were to use to filter's peak transmission instead of this charac- 
teristic value, the tabulated emission-line fiuxes and equivalent widths would all be ~ 12% 
smaller. 

The error bars quoted above for the space density of Lya emitters represent only the 
statistical uncertainty of the fits. They do not include the possible effects of large-scale 
structure within our survey volume. Specifically, if the linear bias factor for LAEs is two (see 



Gawiser et al.ll2007l . for an analysis of the objects' clustering) then the expected fiuctuation 
in the density of Lya emitters measured within a ~ 10^ ^7o~^ Mpc^ volume of space is ~ 30%. 
This value should be combined in quadrature with our formal statistical uncertainty. 

Since Lya galaxies have been observed at a number of redshifts, it is tempting to use our 
data to examine the evolution of the LAE luminosity function. Unfortunately, the samples 
obtained to date are not yet robust enough for this purpose. An example of the problem is 
shown in Figure [121 which compares our cumulative luminosity function (and our Schechter 
fit for a = —1.5) to two measures of Lya galaxies at z = 5. 7. As the figure illu s trates , 



there are large differences between the measurements. If the iMalhotra fc Rhoadd (120041 ) 
luminosity function is correct, then LAEs at z = 3.1 are a factor of ~ 2.5 brighter and/or 
more nume rous than their z = 5.7 counterparts. However, if the z = 5.7 LAE luminosity 



function of IShimasaku et al.l (l2006l ) is correct, then evolution is occurring in the opposite 



direction, i.e., the star- formation rate density is declining with time. Without better data. 
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it is difficult to derive any conclusions about the evolution of these objects. 



Figure [T^ also plots our data again st the predictions of a hierarchical model of galaxy 
formation (ILe Delliou et al.ll2005l . l2006l ). As this comparison demonstrates, our luminosity 
function for z = 3.1 LAEs lies slightly below that generated by theory. This is not sur- 
prising: one of the key parameters of the model, the escape fraction of Lya photons, was 
set using previous estimates of the density of 2; ~ 3 LAEs. Unfortunately, these measure- 
ment s were based on extrem ely small samples of objects, specifically, n i ne z = 3.1 emitters 
from iKudritzki et al.l (120001 ) and ten z = 3.4 LAEs from ICowie fc Hul (119981 ). Since these 
surveys inferred a larger space de nsity of Lyg emi t ters t han measured in this paper, a mis- 
match between our data and the iLe Delliou et al.l (120061 ) models is neither unexpected nor 
significant. 



6. Star Formation Rate Density at 2; ~ 3.1 

Perhaps the most interesting result of our survey comes from a comparison of the galax- 
ies' Lya emission with their i?-band magnitudes. Both quantities measure star formation 
rate: Lya via the combination of Case B recombination theory and the Ha vs. star formation 
relation 

SFR(Lya) = 9.1 x 10"^^ L(Lya) Mq yr^^ (5) 
( Kennicutt 19981 : Hu. Cowie. fc McMahon 1998), and R, via population synthesis models of 



the rest-frame UV (A1570) 

SFR(UV) = 1.4 X 10"^^ Mq yr"^ (6) 



( iKennicuttl Il998l ) . If both of these calibrations hold for our sample of Lya emitters, then a 



plot of the two SFR indicators should scatter about a one-to-one relation. 

Figure [12] displays this plot. In the figure, galaxies with Lya star-formation rates less 
than ~ I.I5M0 yr~^ are excluded by our 1.5 x 10"^'' ergs cm^^ s"^ monochromatic flux 
limit, while objects with large UV star-formation rates, but weak Lya are eliminated by 
our equivalent width criterion. The latter is not a hard limit, since LAE colors range from 
< (-B + V) — i? < 2.5, and it is the B + V continuum that is used to define equivalent 
width. Nevertheless, if we adopt 1.4 as the upper limit on the median color of an Lya 
emitting galaxy (i.e., 1 a above the median {B + V) — R 0.65 color of the population), we 
obtain the dotted line shown in the figure. 



Despite these selection effects, the Lya and UV continuum star-formation rates do 
seem to be correlated. However, there is an offset: the rates inferred from the UV are. 
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on average, about three times higher than those derived from Lya. While the Lya SFR 
measurements are generally less than 10 /^ro"^ Mq yr~^, the rest-frame UV values extend up 
to ~ 50 hm ~'^ Mc^ Thi s discrepancy has previously been seen in a sample of 20 LAEs 

at z = 5.7 (lAjiki et al.ll2003l ). and has two possible explanations. 



The most likely ca use of th e offse t is the galaxies' internal extinction. By studying 
local starburst galaxies, ICalzettil (120011 ) has shown that a system's ionized gas is typically 
attenuated more than its stars. In other words, while optical and IR emission-line ratios 
can usually be reproduced with a simple screen model, the shape of the UV co ntinuum 
requir es that the dust and stars be intermingled. For a self-consistent solution, ICalzetti 
(120011 ) suggests 

EiB - y)3tars = QAAE{B - V)^^, (7) 

If we apply the ICalzettil (120011 ) law to our sample of z = 3.1 Lya emitters, then for the UV 
and Lya star-formation rates to be equal, the extinction within our LAEs must be as shown 
in Figure [TH According to the figure, in most cases it only requires a small amount of dust 
{E{B — y)stars < 0.05) to bring the two indicators into agreement. Figure [H] also suggests 
that internal extinction becomes more important in the brighter gala xies. This is consistent 



with observations of local starburst systems (e.g.. iMeurer et al.lll995l ). and is expected if the 



mass-metallicity relation seen in the local universe carries over to dust content. 

Alternatively, the discrepancy between the Lya and UV continuum star-formation rates 
may simply be due to uncertainties in their estimators. Models which translate UV lumi- 
nosity into star formation rate have almost a factor of two scatter and rely on a number of 
parameters, including the initial mass function and the timescale for star formation. The 
latter is particularly problematic. Lya photons are produced almost exclusively by extremely 
young (< 30 Myr), massive (> IOMq) stars which ionize their surroundings. It therefore 
registers the instantaneous star-formation occurring in the galaxy. Conversely, continuum 
UV emission (at 1570 A) can be produced by populations as old as ~ 1 Gyr; thus, it is 
a time-averaged quantity. If the star-formation rate in our Lya emitters has declined over 
time, then it is possible for UV measure ments to systematically overestimate the present day 
star formation (IGlazebrook et al.lll999l ). 



If we assume that Lya emission is an accurate measure of star-formation, then it is 
possible to integrate the Schechter function to estimate the total contribution of LAEs to the 
star- formation rate density of the z = 3.1 universe. We note that this procedure does carry 
some uncertainty. If we just consider galaxies brighter than our completeness limit (1.5 x 



10 ergs cm ^ s ^ or LLya > 1.3 x 10' 



'^^ hjQ ^ ergs s ^) then the star- formation rate density 



associated with LAEs is ~ 3.6 x 10'^ HjoMq yr"^ Mpc"^, or 1.2 x 10''^ hjo Mq yr"^ Mpc'^ 
if the internal extinction in these objects is E{B — V^)stars ~ 0.05. However, to compute the 



- 18 - 



total star-formation rate density, we need to extrapolate the LAE luminosity function to 
fainter magnitudes, and even 160 objects is not sufficient to define a to better than ~ 25%. 
Consequently, our data admit a range of solutions. 

This is illustrated in Figure [151 which displays SFR likelihoods derived from the prob- 
abilities illustrated in Figure [TTl As the figure shows, the most likely value for the LAE 
star- formation rate density of the 2; = 3.1 universe (uncorrected for internal extinction) is 
6.5 X 10~^ h-jQ Mq yr~^ Mpc~^, while the median value of this quantity (defined as the point 
with equal amounts of probability above and below) is 8.6 x 10~^ /i7o yr~^ Mpc~^. More- 
over, these numbers are likely to be lower limits: if the discrepancy seen in Figure [13] is due 
to internal extinction, then the true SFR density is probably ~ 3.5 times higher. 

The numbers above indicate that at z = 3.1, the star-formation rate density asso- 
ciated with Lya emitters is comparable to that found for Lyman-break galaxies. Before 
correcting for extinction, our number for the LAE star-formation rate density is 8.6 x 
IQ-^ho Mq yr-i Mpc'^ For comparison, the L BG star-formation rate density at z = 3.1 



(before extinction) is ~ O.Olh^o Mq yr'^ Mpc"^ flMadau et al.lll998l : ISteidel et al.lll999l ). It 
is true that internal extinctio n within Lyman-b reak galaxies is typically lar ger than i t is in 
our LAEs, E{B -V) 0.15 Jsteidel et al.lll999h . However, according to the lCalzettil J200lh 
extinction law, the effect of dust on the emission line flux of a galaxy is much greater than 
that on the stellar continuum. Consequently, our dust corrected SFR density for LAEs, 
~ 0.03/i7o Mq yr~^ Mpc~^, is ~ 75% of the LBG value. Of course, given the extrapolations 
and corrections required to make this comparison, this number is highly uncertain. 



Discussion 



The space density of z = 3.1 Lya emitters shown in Figure [TT] translates into a sur- 
face density of 4.6 ± 0.4 arcmin"^ per unit redshift interval above our c omple teness limit. 
This number is similar to that derived by iThommes fc Meisenheimerl (120051 ). under the 
assumption that the LAE phenomenon is associated with the creation of elliptical galax- 
ies and spiral bulges. It is also consis tent with the semi-analytical hierarchical structure 



calculations of iLe Delliou et al. 
ber of z 



(l2005l ). though the latter predict a slightly larger num- 



3 LA 



Le Delliou et al. 



is th an found in this paper. This difference is not significant, since the 



2005h models ha ve be en adjusted t o matc h the previous small- volume Lya 
surveys of iKudritzki et al.l (I2OOOI ) and ICowie &: Hul (119981 ). A ~ 30% re-scaling of the es- 



cape fraction of Lya photons solves the discrepancy, and maintains the match between the 
predictions and the faint-end slope of the galaxy luminosity function. 
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More notable is the excellent agreement between the lLe Delliou et al.l (120061 ) simulations 
and the observed distribution of Lya equivalent widths (Figure [9]). Both are very well- fit via 
an exponential with a large (~ 75 A) scale length. Moreover, the models also predict that 
the scale length observed for a magnitude-limited sample of galaxies (such as that produced 
by the Lyman-break technique) will be much smaller th an that found via an emission-line 
survey. This is consistent with the LBG results found by lShapley et al.l (120031 ). 



Nevertheless, we should emphasize that the LAEs detected in this survey are probably 
not primordial galaxies in their initial stages of star-formation. Very few of the objects have 
the extremely high equivalent widths calculated for stellar populations with top-heavy initial 
mass functions. More importantly, the scatter in the galaxies' mioeo — ''^1570 colors, along 
with the offset between the Lya and UV continuum star-formation rates, suggests that these 
objects possess a non-negligible amount of dust. The existence of this dus t argues against 
the Pop III interpretation of 2; ~ 3 Lya emitters (j Jimenez fc Haimarul2006l ). 



The extremely strong line emission associated with LAEs makes these objects especially 
suitable for probing the evolution of galaxies and structure in the distant universe. The 
space density of z = 3.1 emitters shown in Figure [TT] translates into a surface density of 
4.6 ± 0.4 arcmin"^ per unit redshift interval above our completeness limit. This, coupled 
with our measured luminosity function, implies that in the absence of evolution, there are 
~ 12 LAEs arcmin"^ brighter than 1.5 x 10~^^ ergs cm~^ s~^ in the redshift range 2 < ^ < 4. 



Wide field integral field units, such as those be ing designed for ESO (IHenault et al.ll2004j ) and 
the Hobby-Eberly Telescope (IHill et al.ll2006l ) will therefore be able to find large numbers of 
Lya emitters in a single pointing. Moreover, because the faint-end of the luminosity function 
is steep (a ~ —1-5), the density of LAEs goes hnearly with survey depth. Dropping the flux 
limit by a factor of two (to 7.5 x 10^^^ ergs cm~^ s^^) will roughly double the number of 
LAEs in the sample. 



With an integral-field spectrograph, it is also possible to increase the sample of high- 
redshift galaxies by identifying objects with equivalent widths lower than our detection 
threshold of 80 A (~ 20 A in the LAE rest frame). However, the gain in doing so is likely to 
be small: according to Figure Lya rest-frame equivalent widths e-fold with a scale length 
of ~ 75 A. If thi s law extrapolates to weaker-lined systems, as suggested by the models of 
Le Delliou et al.l (120061 ) . then most Lya emitters are already being detected, and pushing 
the observations to lower eq uivalent widths wil l only increase the number counts by ~ 20%. 
Furthermore, as the data of iHogg et al.l (119981 ) demonstrate, contamination by foreground 
[O II] objects increases rapidly once the equivalent width cutoff drops below ~ 50 A in the 
observers frame (or ~ 12 A in the rest frame of Lya). Unless one can accept a large increase 
in the fraction of contaminants, surveys for high-redshift galaxies need to either stay above 
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this threshold, or extend to the near-IR (to detect B./3 and [O III] A5007 in the interlopers). 
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Table 1. Log of Narrow-band Observations 





Exposure 


Seeing 


Active 


UT Date 


(hr) 


(") 


CCDs 


6 Oct 2002 


2.0 


1.4 


8 


12 Oct 2002 


1.7 


0.9 


8 


4 Jan 2003 


2.0 


1.0 


8 


5 Jan 2003 


3.0 


1.0 


8 


6 Jan 2003 


3.0 


1.1 


8 


29 Nov 2003 


2.5 


1.0 


7 


1 Dec 2003 


1.7 


0.9 


7 


23 Jan 2004 


2.3 


1.3 


7 


24 Jan 2004 


2.0 


0.9 


7 


25 Jan 2004 


2.5 


1.1 


7 


16 Feb 2004 


1.1 


1.0 


7 


17 Feb 2004 


0.8 


1.1 


7 


18 Feb 2004 


1.3 


1.0 


7 


19 Feb 2004 


1.2 


0.9 


7 


20 Feb 2004 


0.9 


1.0 


7 
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Table 2. Candidate Lya Emitters: The Statistically Complete Sample 



ID 


a(2000) 


5(2000) 


Log F5000 


Equivalent Width 


1^ 


03:33:16.86 


-28:01:05.2 


-15.596 


449 


2 


03:33:12.72 


-27:42:47.1 


-15.832 


392 


3^ 


03:33:07.61 


-27:51:27.0 


-15.860 


92 


4 


03:32:18.79 


-27:42:48.3 


-15.888 


251 


5 


03:32:47.51 


-27:58:07.6 


-15.956 


235 


6 


03:32:52.68 


-27:48:09.4 


-15.960 


272 


7 


03:31:44.99 


-27:35:32.9 


-15.972 


248 


8 


03:31:54.89 


-27:51:21.0 


-15.988 


310 


9 


03:31:40.16 


-28:03:07.5 


-16.040 


116 


10 


03:33:22.45 


-27:46:36.9 


-16.080 


143 



^Candidate AGN 

^Foreground AGN 
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Table 3. Candidate Lya Emitters: Objects Beyond the Completeness Limit 



ID 


q;(2000) 


5(2000) 


Log F5000 


Equivalent Width 


163 


03:33:14.82 


-27:44:09.1 


-16.824 


380 


164 


03:32:08.46 


-27:48:43.5 


-16.824 


445 


165 


03:33:26.22 


-27:46:09.0 


-16.824 


468 


166 


03:33:11.73 


-27:46:51.7 


-16.828 


312 


167 


03:31:26.49 


-27:50:34.3 


-16.828 


209 


168 


03:33:05.64 


-27:52:47.2 


-16.828 


284 


169 


03:31:50.46 


-27:41:15.2 


-16.828 


364 


170 


03:33:17.68 


-27:45:44.5 


-16.832 


109 


171 


03:31:48.98 


-27:53:38.7 


-16.832 


322 


172 


03:33:10.77 


-27:52:41.4 


-16.836 


356 
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Tablc 4. Photometric Uncertainties 



Log F5000 o- (mag) Log F5000 cr (mag) 



-15.30 


0.022 


-16.20 


0.065 


-15.40 


0.024 


-16.30 


0.073 


-15.50 


0.026 


-16.40 


0.082 


-15.60 


0.031 


-16.50 


0.104 


-15.70 


0.033 


-16.60 


0.125 


-15.80 


0.038 


-16.70 


0.158 


-15.90 


0.042 


-16.80 


0.204 


-16.00 


0.050 


-16.90 


0.264 


-16.10 


0.058 


-17.00 


0.329 
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Table 5. Schechter Function Parameters 



Parameter Best Solution Marginalized Values 



Log L/L* (ergs s ^) 


42.66 


42.641°:?^ 


a 


-1.36 


-1.49l°:g 


N{> 1.3 X 10^2 /i7o"^ ergs s"^) Mpc-^ 


1.46 X 10-=^ 


lA&tVA X 10-' 


(t)* (Mpc-3) 


1.28 X 10-=^ 
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4600 4800 5000 5200 5400 
X (Angstroms) 

Fig. 1. — The bandpass of our narrow-band A5000 filter, along with those of the B and V 
filters, which are used to define the continuum. A spectrum of a typical z = 3.1 Lya galaxy 
is overlaid for comparison. Our narrow-band filter isolates the emission line of Lya sources 
with 3.09 < z < 3.13. 



I I I I I ■ I I I I I I I I 



18 



19 



I , 
20 



I I I I I I I 

21 22 



I I I I I I 
23 



I I I I 




24 



25 



26 



50QO (AB) 



Fig. 2. — Excess emission in the narrow-band A5000 filter over the continuum for objects in 
our survey field. The abscissa gives the instrumental A5000 magnitude, while the ordinate 
shows the difference between the sources' narrow-band and B + V continuum AB magnitudes. 
Our narrow-band completeness limit of 1.5 x 10^^^ ergs cm~^ s~^ is represented by a vertical 
line; our equivalent width limit of 90 A is shown via the horizontal line. The curve shows 
the expected 1 a errors in the photometry. Candidate emission line galaxies are denoted as 
blue circles; the green dots indicate LAE candidates found by our detection algorithms, but 
rejected upon visual inspection. 
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Fig. 3. — Narrow-band A5000, B + V, and difference images for three candidate emission-line 
galaxies. Each frame is 10" on a side, with north up and east to the left. The objects span a 
range of brightness from logFsooo = —15.60 at the top to logF5ooo = —16.74 at the bottom. 
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Log Flux 



A Log Flux 



Fig. 4. — The left-hand panel shows the transmission curve for our narrow-band A5000 filter 
at the outside ambient temperature and in the converging f/3.2 beam of the 4-m telescope. 
Note that the bandpass is nearly Gaussian in shape. The arrow shows the transmission value 
used to translate AB magnitude into monochromatic flux (see text). The center panel uses 
the transmission ftmction to illustrate how our survey volume changes with emission-line 
sensitivity. The right-hand panel translates the transmission function into the photometric 
convolution kernel that is described in the text. 
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Fig. 5. — The sky coordinates of the 160 candidate z = 3.1 LAEs brighter than our complete- 
ness hmit plotted over our narrow- band 5000 A image. The size of each circle is proportional 
to Lya luminosity, with the largest circle representing 1.25 x 10'^^ /^to"^ ergs s^^. The green 
regions show areas of the chip near bright stars that were excluded from the analysis; the 
large rectangle is the GOODS field. 
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Fig. 6. — The B — R (rest frame mioeo ~ ''^1570) color-magnitude diagram for LAEs. The 
sohd circles repre sent galaxies which have been spectroscopically confirmed as Lya emitters 
( iLira et al.l 120071 ): the cross indicates the lone AGN. The long-dashed line at i? = 24 
represents the typical magnitude limit of LBG spectroscopic surveys; the short-dashed line 
a,t R = 25.5 gives the photometric limit of most LBG observations. Note that the median 
color of our LAEs is quite blue; this is consistent with models for galaxies with recent star 
formation. Note also the large range of colors displayed in the figure. This scatter is greater 
than that expected from the photometric errors, and suggests that the LAE population is 
not homogeneous. 
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Fig. 7. — The U — V versus V — R colors of our z = 3.1 LAEs. The sohd circles show spec- 
troscopically confirmed LAEs, the open circles represent sources observed with insufficient 
signal-to-noise for classification, and the crosses are objects with no spectroscopy. The dots 
are the entire 84,410 object catalog. The polygon is the LBG selection region; the sold curve 
is the track of an LBG templat e spectrum. This tr ack falls inside the selection region in 
the redshift range 2.8 < 2; < 3.4 (jshaplev et allboosl ). Although most LAEs have LBG-like 
col ors, their R > 25.5 m agnitu des exclude them from the "spectroscopic" samples studied 
by lSteidel et al. Jl996al jbl. [2003l l The contribution of each LAE's emission line to its l^-band 
flux has been subtracted to yield Vcorr- 
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Fig. 8. — The Rab (rest frame 1570 A) luminosity function of our z = 3.1 Lya emitters (solid 
circles), compared to the res t-frame 1700 A lum inosity function of z = 3.04 Lyman-break 
galaxies (open circles) from (jSteidel et al.lll999l ). The flattening of our luminosity function 
at i? > 26.5 is due to selection: at these magnitudes, only the strongest line emitters make 
it into our sample. In the magnitude range R < 25.5, z = 3.1 Lya emitters are ~ 3 times 
rarer than Lyman-break galaxies. 
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Fig. 9. — The top panel shows the observed distribution of equivalent widths for all the Lya 
emission-line galaxies in our sample. The dotted line shows the apparent best-fit exponential 
for the distribution; the solid curve shows the exponential after correcting for the effects of 
photometric error and our filter's non-square transmission curve. The lower two panels divide 
the sample in half, and demonstrate that the exponential law does not change much with 
galaxy luminosity. The vertical dashed line shows the maximum equivalent width expected 
for populations with normal initial mass functions. 
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Fig. 10. — The number density of z = 3.1 Lya galaxies with observers' frame equivalent 
widths greater than 90 A binned into 0.2 mag intervals. The points give the density of 
objects under the assumption that our filter's FWHM defines the survey volume; the open 
cir cles repres e nt dat a beyond our completeness limit. The solid curve shows our input best- 
fit ISchechterl (119761 ) luminosity function, while the dashed line illustrates the shape and 
normalization of this function after correcting for the effects of photometric error and our 
filter's non-square transmission curve. 
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Fig. 11. — Maximum likelihood confidence contours for our ISchechterl (119761 ) function fit 
to the observed distribution of Lya fluxes. The three parameters in the analysis are the 
faint-end slope (a), the bright-end cutoff (logL*) and the space density of galaxies with 
observed monochromatic fluxes greater than 1.5 x 10^^^ ergs cm~^ s~^, i.e., Li^ya > 1-3 x 
10^^ hjQ~'^ ergs s~^. The contours of probability are drawn at 1 a intervals. 
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Fig. 12. — The cumulative Lya luminosity function inferred from our survey of Lya emitters 
with rest- f rame equivalent widths greater than 22 A. The solid blue line shows our best fit 

-1.49), the green line is the a = —1.5 luminosity function 
Malhotra fc Rhoadsl (120041) for LAEs at z = 5.7 , and the red line is the Schechter 



found by 
fit for z : 



Schechter (|l976l) function (a 



5.7 ernitters found by IShimasaku et al.l (120061 ). The dashed line is Model A by 



Le Delliou et al.l (120061 ). For purposes of this figure, our data have bee n artificially normal- 



i zed t o mat ch our best-fit function . The large difference between the iMalhotra fc Rhoads 
( I2OO4J ) and IShimasaku et al.l (120061 ) fits makes it impossible to study the evolution of the 
LAE population at this time. 
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Fig. 13. — A comparison of the star formation rates derived from Lya emission (under 
Case B recombination) and t he UV continuum at 1570 A. The sohd dots represent spectro- 
scopically confirmed objects (ILira et al.l 120071 ). The diagonal dashed hne shows where the 
two measurements are equal, while the solid line illustrates where the UV continuum star- 
formation rate is three times the Lya rate. Our flux limit is shown via the vertical dashed 
line; the approximate location of our equivalent width threshold is shown via the dotted line 
(see text). Note that, although the two indicators are correlated, the Lya-inferred rates are 
~ 3 times less than those derived from the UV continuum. 
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Fig. 14. — Estimates of the internal extinction within our Lya emitters, fo rmed using th e 
assumption that the galaxies' ionized gas is attenuated more than its stars (jCalzettill200ll ). 
The solid dots represent spectroscopically confirmed Lya emitters; the dashed line shows the 
monochromatic flux limit of our survey. Note that very little internal extinction is needed 
to bring the UV and Lya star- formation rates into agreement: even in the bright {R > 24.5) 
galaxies, the internal extinction is never more than E{B — V^)stars ~ 0.1. 
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Fig. 15. — The results of our maximum likelihood analysis for the contribution of Lya 
emitters to the star formation rate density of the universe. The abscissa is the star formation 
rate density derived from the observed luminosity of the Lya emission line; the ordinate is 
the relative probability of a solution. The dashed line shows the observed star formation rate 
density associated with galaxies above our completeness limit, i.e., without any extrapolation 
of the galaxy luminosity function. No correction for internal extinction has been applied. 
The figure implies that the amount of star formation taking place in galaxies with strong 
Lya emission is comparable to that in Ly man-break galaxies. 



